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Abstract 

The work deals with one of the topics of collective motion. In the frame 
of Nuclear Fluid Dynamics, a model which portrays the nuclear matter as a 
quantum elastic body, the torus-like motions and their associated energies 
■ are computed using the thirteen moment approximation. Such excitations 

correspond to the Hill vortex known from classical Hydrodynamics. There 
are also calculated the nonvanishing contributions of transverse electric form 
factors and differential cross sections in the electroexcitation of these col- 
lective modes, which are purely toroidal. The spin-dependent collective ex- 
citations, with toroidal electromagnetic structure, are studied by means of 
the Generalized Goldhaber- Teller model, with emphasize on the 1~ spin-flip 
mode and its excitation in spherical nuclei by inelastic electron scattering. 
We discuss the importance of toroidal contributions in the inclusive electron 
scattering (e, e') and exclusive coincidence electron scattering (e, e'7). In or- 
der to extract the toroidal multipole, we use the backscattering angles in the 
first mentioned reaction, and the separation of the longitudinal/transverse 
interference in the second case. The introduction of a quantity which ac- 
counts for the deviations from the Siegert theorem, shows the importance 
of toroidal quadrupole transitions at high-momentum transfer. Another im- 
portant result concerns the dependence of the intensity of toroidal effects on 
the nuclear vorticity. 



1 Introduction 

The toroidal multipole moments are a distinct family of electromagnetic moments, 
which occur in a special parametrization of charge and curent densities Q, [@. 

Long time ago, Zeldovich showed that a 1/2 spin particle may interact with an 
electromagnetic field not only by means of a dipole electric interaction d(cr-E) that 
simultaneously violates the space and time inversion, but also by a different type of 
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interaction, a{u ■ J), which violates the parity but not the time reversal (Figure 1). 
This new electromagnetic characteristic was named anapole and it is the first 
member of the class of toroidal moments, i.e. the static toroidal dipole moment. 
Inside nuclei, the anapole moment arises as a consequence of P-noninvariant Q 
nuclear forces. In classical electrodynamics a very simple example of a toroidal 
dipole is given by a solenoid folded onto a torus. The radiation resistence of the 
corresponding antenna is proportional to (-R^t/A 3 ) 2 , where and are the 
torus small and large radii ||. Therefore the toroidal antenna depends on the ratio 
between the geometrical size of the source d and the radiation wavelength A like 
the charge octupole and magnetic quadrupole. Such solenoids have non-vanishing 
magnetic potential in those regions of the space where the toroidal dipole moment 



The work that we present below discusses the problem of those nuclear collec- 
tive motions associated to purely or partially toroidal electromagnetic transitions. 
There will be also studied the properties of excitations associated to vortical nu- 
clear currents of orbital isoscalar nature {dipole torus mode : DTM) and of isovector 
spin dependent nature {spin - isospin mode : s-is) along with the investigation of 
the possibility to detect them in electron inelastic scattering on nuclei. 

The electroexcitation of collective rotational and vibrational motions with the 
account of toroidal quadrupole transitions will be also investigated. 



2 The Dipole Torus Mode in Nuclear Fluid Dy- 
namics 



The Nuclear Fluid Dynamics (NFD) allows the description of isoscalar giant reso- 
nances as in phase harmonic vibrations of proton and neutron fluids, the restoring 
force being a consequence of the elasticity of nuclear matter. The nucleus is por- 
trayed as a Fermi fluid which as a result of an external perturbation may undergo 
longitudinal oscillations (irotational), as is the case of the liquid drop, but also 
transversal oscillations. Consequently an isoscalar giant resonance will be a de- 
scribed by a small amplitude collective oscillation with multipolarity A, in which 
protons and neutrons are performing a divergenceless irotational or vortical dis- 
placement [jj. 

Basically, the NFD equations can be deduced by taking the classical limit of 
time dependent Hartree-Fock equations (TDHF) 



which is nothing else than the collisionless Boltzmann equation from statistical 
physics of transport phenomena : 



^4 = [H, p] 
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and represents the equation of motion for the Wigner transform of the density 
matrix p 

f(r,p,t) = Jdsp(r + ^r-^t)e~^ s (2.3) 

In the argument of / the position vector r of a particle is determined with a 
precision larger than the wavelength of the perturbation in such a way that the 
Heisenberg principle is not violated. 

In order to reduce the complexity of the non-linear partial derivatives of equa- 
tion (2.2), one performs a transformation to coupled equations for the macroscop- 
ical variables. These physical quantities are introduced as p-order moments ( = 
for the density p, — 1 for the three components of the mean velocity Ui(r,t), = 2 
for the nine components of the strain tensor Pij(r,t) ) of the distribution function 
with respect to the momentum as follows : 

p(r,t) = mjdpf (r,p,t) (2.4) 

p (r, t)ui(r, t) = J dp f (r, p, t) p< (2.5) 

Pij(r,t) = — ( dp f(r,p,t)(pi - mui){pj - muj) (2.6) 

By integrating the equation of Boltzmann over the momentum p, with weight 1, 
Pj/m, piPj/m 2 , imposing the incompresibility condition for the nuclear matter, i.e. 
P = Po =constant and assuming for the strain tensor the ansatz 

P ij = Po+Pi j (2.7) 
one get the linearized NFD equations (the thirteen moments approximation) 

du k 



dx k 



(2.8) 



It is convenient to write the displacement dxi of the nuclear matter at a certain 
point r inside the nucleus as 

dxi = a$(r) (2.11) 



or, alternatively by the mean velocity 



«i = a?(r)^ (2-12) 

where a^(r) is the vector field of instantaneous displacements in the Fermi con- 
tinuum, and a\ is the time-dependent amplitude of harmonic oscillations (a\ ~ 
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smu)\t) associated to the resonant phenomena that we study in this paper. If we 
differentiate (2.9) with respect to time, (2.10) with respect to position x%. and we 
use equation (2.8), we obtain 

d 2 Ui d 2 u t 

P ^ = Po M (2 ' 13) 

Afterwards, employing the ansatz (2.12) we arrive at the Helmholz equation (see 
(2.1)) for stationary spherical waves 

^ + k 2 a x = (2.14) 
oxj 



where k = J puj 2 /Pq is the wave number. Equation (2.14) admits three independent 
solutions || 

a.i=N?Vjx(k l r)Y Xli {0,<f>) (2.15) 

a x = N X V x rj x (k t r)Y XfM (9, 0) (2.16) 

a x = N X V x V x rj x (k t r)Y Xfi (9, 0) (2.17) 

in a frame with fixed axis. The longitudinal and poloidal solutions describe com- 
pressional and transversal oscillations of the elastic nuclear globe, being responsible 
for the electric-like resonances with parity n = (— ) A , whereas the torsional solution 
describes magnetic-like resonances with parity n = (— ) A+1 . In the longwavelength 
limit (kr < 1), the poloidal vector field becomes proportional to the longitudinal 
one 

a x (r) =iV p A Vx Vxrr A y AO (#,0) = A^A + l)Vr A r A0 (M) = (A + l)a A (r) (2.18) 
and the torsional field is merely 

a x (r) = iV A V x rr x Y X0 (6,<P) (2.19) 

It is important to substantiate that in the above mentioned limit, the poloidal 
solution is simultaneously irotational and divergenceless, i.e. V-a A = V xa A = 0, 
whereas the torsional solution is purely solenoidal, V-a A =0 but V x a x ~ a x ^ 0. 

In the dipole case (A = 1), the solution (2.18) corresponds to the displacement 
as a whole of the nucleus, without the change of the internal state. The equality 
(2.18) is obtained from (2.15) and (2.17) by keeping the first term in the asymptotic 
expansion of the spherical Bessel function 

^> ^ f 1 - 2(2173) + -") (2 ' 2 °) 

Consequently, in order to investigate the dipole response of an incompressible 
elastic globe one need to go beyond the limit imposed by the longwavelength ap- 
proximation and to introduce the high-order terms in the expansion (2.20). Then, 
the poloidal solution becomes || 

a l p = N l p V x V x rr 3 Y 10 {6, <p) (2.21) 
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In order to establish an explicite expression for the displacements field, corrected in 
such a way to take into account the center of mass motion we impose the condition 



f drpal 

SRcm. = , , P = 2.22 
Jdrp 

This procedure leads to the following expression for the displacements field 
a\ = NpV x V x rr(r 2 - R 2 )Y l0 {6, 0) 

= ^=N l p [V2r 2 Y° 12 (0, <P) + (5r 2 - 3i? 2 ) Y%(6, 0)] (2.23) 
We can rewrite eq.(2.23) on spherical components 

a 1 p )r = \^-N 1 p (r 2 -R 2 )cose (2.24) 



7T 



■pi 




3 M(2r 2 -i? 2 )sin^ (2.25) 



7T 



(oJ), = (2.26) 

The dipole poloidal displacement field, or the Dipole Torus Mode (MDT) coincide 
with that for the Hill vortex known from Hydrodynamics ||10|| . The Stokes current 
function corresponding to this vortical flow is given by 

^(r, 0) = Np (r 2 - R 2 )r 2 sin 2 (2.27) 

The contour lines given by eq.(2.27) are ploted in Fig. 2. The critical or stagnation 
points are fixed by the conditions a r = and ag = 0, i.e. r c = R/V2 §i 9 C = ±vr. 
The geometric locus of these points is represented by a ring in the equatorial plane 
of the sphere. The vortical flow of the fluid takes place around this ring. The 
stream lines rotated around the globe axis generate tori. Such a vortical flow is 
known in Classical Hydrodynamics under the name of Hill ring vortex, contrary to 
the linear vortex where the critical points are located on the symmetry axis of the 
spheroid. We called the collective excitation corresponding to the Hill vortex of 
the nucleus dipole torus mode ( DTM ) ||. A similar kind of collective motion have 
been studied in [|llj and |12| . In the equatorial plane, i.e. the plane containing the 
vortical ring ( 6 — | ) 



C = (V x u p ) e=K = h\ -Nla(t)re^ (2.28) 
2 y 7r 



Thus, the vorticity depends on the radial coordinate. The energy is given by 

E(\- or ) = \ —hcu F w 2hu (2.29a) 
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Therefore, DTM may be interpreted as a dipole transversal isoscalar resonance of 
2%uj type. Using realistic parameters for the Fermi distribution we get ||13||, 



E(\- or ) = 93.72A- 1/3 MeV (2.296) 

i.e. the predicted mode is most probably located between the giant isovector 
resonance and the giant isoscalar octupole resonance. 

We would like to complete the analisys on toroidal dipole by presenting the 
calculation of the transverse electric form factor and of the current and transition 
vorticity densities using the Born's plane wave approximation. The knowledge of 
these quantities is important because the form facor may be directly measured in 
electron inelastic scattering processes |14| . 



In the case of the DTM, the current density associated to the transition is given 
in the fluid-dynamic representation by 

J tor = n e u p = n e a l p {r)a(t) (2.30) 

where n e = eZ/Ano and hq = 3A/4nR 3 is the particle density ; a(t) = aoucosuit 

and a = (dfe)^ ( 5 i = ^( N p) 2MR ^ C ^ = ^( N l? Mv F R2 )- We introduce the 
electric transverse form factor 

Att 

Wt\k)\ 2 = -{\K{k,t)\\ (2..31) 

where by < ... > t we understand the time averaging. Normalizing to c 2 one obtains 
the dimensionless quantity 



c 2 



In Fig. 3 we represented \Ff(k)\ 2 for three spherical nuclei : 40 Ca, 90 Zr and 208 Pb. 
This plot exquibite an enhancement of the DTM in heavy nuclei compared to light 
nuclei. Moreover, we notice that the first diffraction maxima of the form factor is 
shifted towards small momentum transfer when we pass to nuclei with Z and A 
large. 

A simple calculation shows that the longitudinal and magnetic multipoles 
vanish for the current density (2.30). Thus, the transverse electric form factor has 
the only non-vanishing contribution in the excitation with electromagnetic probes 
(photons,electrons) of DTM! The Coulomb multipoles does not participate in the 
excitation due to the fact that the charge transition density p^ or vanishes as a a 
consequence of the incompresibility condition (2.8) imposed to the Fermi globe. 
In other words, since the excited mode is purely rotational, the longitudinal form 
factor vanishes. Although the magnetic multipole should be associated to the 
excitations of rotational motions, these are of oposite parity to DTM, which has a 
natural parity like the electric excitations. 

The electroexcitation differential cross-section of DTM looks as follows 

= ^a Mott f-l ( + tan 2 ~) || ff(q) || + )| 2 (2.33) 



1 tor 



2q 2 2 
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The dependence of (2.33) on the scattering angle is given in Fig.4. We should 
notice that since the longitudinal part does not contribute to the electroexcitation 
of DTM, it does not appear necessary to consider the backscattering case. However, 
in the electroexcitation process there will occur also modes which contribute to the 
longitudinal part and consequently in order to separate their influence it is better 
to chose the case with scattering angle 9 = 180°. 

An important characteristic in the study of collective oscillations by inelastic 
electron scattering is given by the transition current density , a quantity suscep- 
tible to experimental determination. Basically we are interested in the multipole 
component J~\\+i(r) which can be expressed as inverse Fourier transform of the 
transverse electric form factor [|15[| : 



f°° 

J 12 (r) = --= / Ff(q)j 2 (qr)q 2 dq (2.34) 
which integrated gives 



Jl2(r) = -*fciF- 0<r<R 

= 7=-^, r = R 

= 0, r>R 

This radial function is ploted in fig. 5. 

Another quantity of interest in (e, e') inelastic processes is the vorticity 
It determines the nuclear current properties unconstrained by the charge-current 



conservation law. The vorticity transition density is [16 



oo 11 (r) = Vs(J- + ^)j 12 (r) (2.35) 



This equation tells us that inside the nucleus the vorticity density varies linearly 
with the radius in the same manner as the vorticity vector (2.28). This radial 
function is represented in Fig. 6. 

Introducing the multipolar parametrization of Dubovik and Cheshkov the 
transverse electric multipole corresponding to DTM reads 

= <tft (2-36) 

if we take into account that the Coulomb multipole Qio — 0. For the toroidal 
dipole moment associated to the transition + — > \^ or the following proportionality 
relation is available 

(7I) t ~ aZ (2.37) 

A similar result was known for the electron transition Is 1/2 — > 2pi/ 2 in Hydrogen- 
like atoms JTJJ. This is the reason why the electromagnetic effects of toroidal nature 
are so weak for small Z. But whereas in the above mentioned atomic transition 
the toroidal moment enters as a small correction, in the DTM transition it is the 
principal electromagnetic characteristic of the nuclear response. 
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3 The spin-flip resonance 



If we consider the excitation of nuclei with + ground state, then the Wigner 



supermultiplet theory [IS leads to the classification of giant dipole resonances 
given in Table 1. In the frame of generalized Goldhaber- Teller mode for isobaric 
nuclei (N = Z), having ground state with J = + , T = 0, the current density 
associated to the GT mode is 



J{r) = ^QnPoir) 



(3.1; 



and the spherical components of the magnetization density, for the s-is and sw 
modes, are given by 



{ft 



{ft 



h 9p-9n 

4mc 2 

h g P + g n 



(3.2) 



(3.3) 



4mc 2 

where Po{r) is the ground state density, and q n is the relative coordinate. The 
contribution of the sw state (3.3) may be neglected due to the smallness of the 
factor [(g p + g n )/ (g p - g n )} 2 . 
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The calculation of electromagnetic multipoles shows that whereas the GT mode 
is mainly longitudinal the s-is one is purely toroidal like DTM. It is worthwhile to 
make a comparison between the electroexcitation differential cross-sections of the 
isovector resonances 1~ GT and 1~ s-is. 




CMottb' 



l-GT 



1 s— is 



2A 



q 2 V L {6) + 2 



F 2 (q) ( q 2 kc\ fg p -g n 



OMottb \^— 2 j 



V T (9) 



V T {9) 



(3.4) 
(3.5) 



where b = ^Jmuj/h is the characteristic length of the harmonic oscillator with 
frequency u, and the charge density form factor in the ground state is 

F(q) = I dre^ r p (r) (3.6) 
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In Fig. 7 we have ploted the differential cross-sections of the two resonances, GT 
and s-is, corresponding to the electroexcitation of the 12 C and 16 nuclei. Notice 
that for small scattering angles, the differential cross-sections of the GT resonances 
are much larger than those for the s-is resonances, even at large momentum trans- 
fer. However, for backscattering angles, the differential cross-sections of the s-is 
modes become important and exceed those of GT resonances at momentum trans- 
fer q > 0.5fm _1 . 



4 Toroidal quadrupole transitions in the Riemann 
Rotational Model 

The dynamic character of the nuclear rotational motion is one of the basic prob- 
lems, unsolved yet, in nuclear structure theory. Many attempts have been made 
up to the present time to clarify whether the nuclear matter can be portrayed as 
a quantum fluid that could support only irotational flows (IF) or it is a quantum 
rotor which gives rise to a rigid rotation (RR) of the whole nucleus. In a naive 
image irotational flow may be viewed as a deformation which propagates on the 
surface of the nucleus, along with a corresponding motion of the intrinsic structure 
with small angular momentum (Fig.8c). At the other extreme lays the rigid body 
rotation (Fig. 8a). 

The calculation of inertia moments of these two different types of flow un- 
derestimates the predictions of IF and overestimates those of RR, in comparison 
with the experimental values. Thence one must consider a model which takes into 
account the existence of currents with intermediate values between IF and RR. 



The Riemann rotational model [19] is a simple generalization of the Bohr- 



Mottelson model, with the current ranging between the limits mentioned above. 
Moreover the associated velocity field is supposed to depend linearly on position. 
A Riemann rotator is an elipsoid whose principal axes have lengths 01,02,03 and 
in stationary conditions are at rest with respect to a frame rotating with constant 
angular velocity u). In this rotating frame internal motions with vorticity C, = 
Vxu occur. These two vectors are parallel and oriented along one of the ellipsoid 
principal axes. Thence the velocity field measured by an observatory at rest with 
respect to the rotating frame is given by 



ui = 2 Ql 2 C3^2 + o" 1 2 (2^3 (4.1a) 



a a 

'"2 = — 2 , 2 + 2 , 2 ( 4 - 1& ) 
a-2 + 03 a 2 + a{ 

0? a? 

u 3 = — 2 3 2 ( 2 :ri + 2 , 3 2 (ix 2 (4.1c) 
a{ + a\ a\ + a 2 

whereas the velocity field defined in a space fixed inertial frame and projected on 
the rotating frame 

U = u + u>xx (4.2) 
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A Riemann sequence is characterized by a parameter 

/ = - - 2 (4.3) 

on 

When / = 0, the ellipsoid is rigidly rotating, and when / = —2 the flow is 
irotational. Defining the rigidity parameter through the relation 

r = l + f - (4.4) 

the RR case is reproduced for r = 1 and the IF case for r = 0. Considering that 
£ and to are parallel with the principal axis x%, and that a\ > a 2 , the velocity field 
(4.2) reads 

a 2 - a 2 

U(r) = (1 - r)-| |cjV(xiX 2 ) + rcj x r (4.5) 

The above equation states that the velocity field U(r) is a convex combination of 

2 2 

rigid U rr = uj x r and irotational L/jf = °a + °I ^V(xiX 2 ) contributions. A similar 
relation is valid for the inertia moment 1 r . A particular example of Riemann flow 
r G (0, 1) is represented in Fig. 8(b). 

In order to check the predictions given by the Riemann model a direct deter- 
mination of the nuclear current is required. The electron-nucleus scattering is an 
usefull tool which allows the measurement of the electromagnetic charge and cur- 
rent densities inside the ground-state band PU| . We showed in a previous section 



that Coulomb multipoles are associated to the charge distribution of the nucleus 
and thus, in order to obtain the quantities depending on the nuclear current we 
need to determine the transversal part of the cross-section, i.e. the electric and 
magnetic multipoles. 

Since the task of this work is to substantiate the toroidal multipoles which are 
active in such processes we will foccus on the study of electric transverse multipoles 
arising in the Dubovik and Cheshkov multipolar parametrization. 

Let us consider an even-even nucleus whose surface oscillates harmonically and 
simultaneously undergo a rigid rotation around an axis perpendicularly on its 
symmetry axis. Expressing the velocities Urr and XJ if as one and two-rank 
spherical tensors, we readily obtain the spherical component \i of the total velocity. 

U lfl (r) = (1 - r) [V 2 ® n] lfi + r[V 1 ® (4.6) 

where the spherical components of the tensors being coupled to = ^/4^rYi M are 
given by 

V -L-RR 

with /J, — ±1, ujfj, — — and T IF and T RR are the inertia moments of IF and RR 
models. 

The current density reads 

J(r) = ff(r)U(r) (4.8) 
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where 



^ = £ p£(r)lio(0, 

L>2 



(4.9) 



is the proton charge density expanded in even multipolar components ( L = 2,4, ... 
) of an axially symmetric nucleus. 

We shall express the electric transverse (2.4) and longitudinal (2.3) multipoles 



as follows [21 
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v 2A + 1 VL 
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A' 






0J 
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1 





. A A 



(4.11) 

where ' ' '] and < ' ' > are 3j and 6j coefficients [E2|. In the RR case we 



\jl —\l —Mi 



take k = 1 and an axial symmetric quadrupole static deformation in (4.9) : /3 
fa + 0, 7 = 



4?ri2§ 

Thus the charge quadrupole components of the density reads 



ff 2 (r)= dnY* (e,<j>)if(r) 



3eZ{3 



5(Ro - r) 



and we obtain in place of(4.10) and (4.11) 



T°l(q,RR) = -W* 



7T- 



40 R 



ji(qRo) - -jsiqRo) 



L 2tM (q, RR) = _Ze ^ 7r "^/f^ b'i(^o) + j 3 (qRo)]^^ 



(4.12) 
(4.13) 

(4.14) 
(4.15) 



For the irotational case = 2 we consider the deformed charge density distribution 
with monopole component 



Po( r ) 



3eZ 



(4.16) 



such that the multipoles (4.10) and (4.11) may be rewriten 



T*l(q,IF) = -ZeVS 
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(4.17) 



L 2M (g,IF) = L 2M (g,RR) (4.18) 

Using the convexity property of the velocity field (4.6) the electromagnetic multi- 
poles for a certain value of the rigidity parameter are 

f^(q,r) = -ZeV3n^^ [ji{qRo) - (l - jjr) j 3 (qR )\ ^ (4.19) 




2g 



L 2At (g,r) = -ZeV2n^—J-^ [ 3l (qR ) +j 3 (qR )}^ (4.20) 



where Qq = ^J^Rq[3 is the static quadrupole moment, and Rq = TqA 1 / 3 . The 
above formulas allow a first interesting remark : The longitudinal multipole does 
not depend on the rigidity parameter and is proportional to the transverse electric 
multipole in the IF limit when r = 




L 2 »(q,r) = \l-T^(q,r = 0) (4.21) 

Therefore the longitudinal multipoles are insensitive to rotational components of 
the velocity field, their values being constant for any value of r. 

Next we will foccus on the behaviour of electromagnetic multipoles at small 
momentum transfer. 

In 1937, Siegert showed that using the charge- current conservation law in the 
longwavelength approximation it is possible to replace the charge density operator 
in the expression of the electric transverse operator with the charge density rate. 
Quantitatively this theorem may be expressed as follows 

"-tell 




n l (q^0) = \ —L x (q^0) (4.22) 



Thus, independently of the model used for the nuclear current the transverse elec- 
tric multipole is proportional to the longitudinal one in the longwavelength approx- 
imation. The approximate equation (4.21) is similar to the exact equation (4.22). 
In other words the Siegert theorem, quantitatively expressed by (4.22) is valid in 
any order of q for the irotational Riemann sequence r = 0. This means that the 
reactions performed at low-g are not able to provide informations on the vorticity. 
The Riemann rotator behaves like an irotational liquid drop at low momentum 
transfer. From the view point of the multipolar parametrization adopted in this 
paper this fact may be justified invoking the following argument : in the low-q 
limit 

(// || L\(q -> 0) || h) = -^(I f || 4(<? - 0) || h) (4.23) 

where Q 2fl = J drr x Y\(9, <f))p(r,t) is the charge quadrupole moment. Using the 
charge- current conservation law, the time derivative of the charge quadrupole op- 
erator may be writen as follows 



Q 2il = / drr x Y^{6, 0) V • J(r, t) (4.24) 
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The presence of the gradient operator in the above equation ensures the cancel- 
lation of rotational (vortical) components of the nuclear current. Thus the time 
derivative of the charge quadrupole moment describes the curless quadrupole flows, 
i.e. V x J = 0, like the well known f3 and 7 vibrations. This is the reason why the 
nuclear response is vibrating-like, without shear components for small momentum 
transfer. Consequently, in order to get informations on the rotational currents in- 
side the nucleus, one needs to investigate the electromagnetic structures negleted 
by applying the Siegert theorem. To go beyond the limitations of this theorem 
one needs to increase the energy transferred in the scattering reaction. The higher 
order terms in the q 2 expansion of the electric transverse multipole are free of the 
constraint dictated by the continuity law and therefore they are likely to provide 
data on the vortical components of the current V x J / 0. Applying the Dubovik 
- Cheshkov parametrization [2] and defining a quantity which takes into account 
the deviation from Siegert's theorem 



V2(q) = q 2 



Tl or (q) || h) 



(If || Tf (0) || Jj) 



(4.25) 



we are in the position to describe the importance of the toroidal multipoles [23] 
In the IF and RR cases, the function r] 2 looks 



%(?,IF) = 
772(9, RR) 



15 



15 



(3 — (gi? ) 2 ) sing.Ro — 3qR cos gi? c 



- 1 



{qRlf 



{{qRi 



) 2 - 2) sin qR + ( 2 - ) qR cos qR 



(4.26) 
- 1 
(4.27) 

whereas in the intermediate case, this function will be expressed as a convex com- 
bination of IF and RR contributions 



772(9, r) = g 2 



(If\\T^(q,r) 
{If II Ti l (0) || Jj) 



1 - r)rj 2 (q, r = 0) + rrj 2 (q, r = 1) (4.28) 



In Fig. 9 we draw this quantity versus momentum transfer q. Notice that the RR 
model presents a stronger deviation from the Siegert approximation than the IF 
model. In the hexadecupole case ( A = 4 ) the deviation effect is smaller than in 
the quadrupoe case ( A = 2 ) for both models. 

Another interesting quantity is the real electric transverse form factor which in 
the lowest order of q 2 is 



mq) 



(I f I) T?(q) [I h) 
(If || Tjf(0) || Jj) 



1 



q'% 

3 r Q 2 



(4.29) 



where Q 2 = 3e 2 Z Rq(3 j 'An is the transition charge quadrupole moment, and T 2 is 
the transition toroidal quadrupole moment. This last equation allows the mea- 
surement of the transition toroidal quadrupole moment T 2 similar to that of the 
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charge (magnetic) mean square radius < r 2 >c(mag) [24], which consists in the com- 
putation of Tfi^q). For an arbitrary value of the rigidity parameter the toroidal 
quadrupole moment reads [[25 



%S 



3e 2 Z3 + 2r 
56tt Z r KoP 



(4.30) 



We ploted in Fig. 10 T 2 versus the rigidity for 152 Sm and 166 Er. An important 
conclusion that we draw from this figure is is that T 2 increases sharply for values 
of r close to 1, being two orders of magnitude larger in the RR case than in the 
IF case. We conclude that in those nuclei where the vortical components of the 
nuclear current are stronger the toroidal quadrupole transitions are intensified with 
respect to the spherical vibrating nuclei. 

Taking into account the definition of the toroidal multipoles, their measure- 
ment is equivalent with the measurement of electric transverse multipole at high 
momentum transfer followed by the removal of the Siegert limit. 

In the excitations of ground state band of an even-even nucleus there will be 
involved the electric multipoles with A = 2,4, ... and and the magnetic multipoles 
with A = 1, 3, .... Further we will not discuss the magnetic multipoles and we will 
foccus on the longitudinal and electric transverse parts of the differential cross- 
sections. Since our purpose is to shed light onto the nature of toroidal transitions 
and their connection with nuclear vortical currents we will chose a convenient 
method to separate the dominant longitudinal components in the differential cross 
section and eventually to obtain the transverse multipoles. Separating the trans- 
verse multipoles by 180° scattering is a seducing method since at this angle these 
multipoles dominate the differential cross section. Another method of separation 
consisting in the scattering of polarized electrons will be discussed bellow. 

For the + — > A + excitation we have the following expression for the differential 
cross-section: 



ott I y/x I / \ + II »> f„\ II n+\|2 i / V i + „„2 ° \ |/\ || r^el 



^) (e,e') f r 





Mx(q) || + >| 2 + + tan 2 - |(A || T A e( (g) || 0" 



(4-31) 

Using the Dubovik - Cheshkov decomposition in the above formula and neglecting 
the Coulomb multipoles and the low-g limit of electric transverse multipole one 

(g (v ,r^ 2 (S +tm i) i<AiiTrwiio+>|2 {m 

This approximation is equivalent to the neglecting of transition charge moments 
Qx mean 2n-power of charge distribution radii r 2n . 

The differential cross sections (4.31) and (4.32) are ploted in both cases, for 
RR and IF, in order to compare the exact formula and the approximate one in 
backscattering processes. The cross section of the quadrupole transition induced 
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by the scattered electron is represented in Fig. 11 for the 152 Sm nucleus. The 
hexadecupolar transition is considered in Fig. 12. From the study of these two 
graphics it becomes obvious that the above mentioned approximation is acceptable 
in the RR model for the whole range of momentum transfer for both considered 
transitions. The main difference consists in the location of diffraction minima. 
These differences are not important because in a phase-shift analysis the curve 
will be smothed in the neighborhood of minimas. 

In the IF case, the exact differential cross section (4.31) matches quite well the 
approximate one for momentum transfer q <400 MeV/c when A = 2 and q <250 
MeV/c when A = 4 for the case 166 Er. The reason of these discrepancies in the 
IF model between the exact and approximate curves is that at high momentum 
transfer the mean 2n-power charge distribution radius is drastically enhanced. 

5 Separation of toroidal multipoles in electron 
coincidence processes 

As we saw in the preceding sections of this work a method to separate the electric 
transverse multipoles is given by the Rosenbluth decomposition of the (e, e') differ- 
ential cross section at 180° scattering angle. A more recent method which allows 



this separation is based on the (e, e'j) longitudinal/transversal interferences |26 
The differential cross section of the process is given by 

/ d 2 a \ h ' a 1 + + / r A+-(A- 2 )+\ 

where Sq +0+ is the differential cross section of the corresponding (e, e') process, 
given by (4.31), r^ + ^( A - 2 ) + is the photodisintegration width of the transition A + — > 
(A — 2) + and T^ tal is the total decay width for the state A + . The ratio of these two 
widths is close to unity. The differential cross section contains as labels the incident 
electron helicity, h — ±1, and the circular polarization of the photon detected in 
coincidence, a = ±1. The angular distribution functions are normalized to unity ( 
J cLVL^Wy, = 1 ), and satisfy the integral condition ( / dVL^W^ = ). Their explicite 
expressions are 

W s (9„ <f> 7 ) = i- (1 + A E (# 7 , 7 )) (5.2) 

W^ 7 ,0 7 ) = ^4 A (fl 7 ,0 7 ) (5.3) 

In the electroexcitation of the rotational g.s. band of a deformed even-even 
nucleus ( + , 2 + , 4 + ,...), the longitudinal multipoles ( C2, C4,... ) are larger 
than the electric transvere multipoles ( E2, E4,... ). The even-even nuclei are 
not good candidates to become polarized targets, because they have spin and 
parity + in the ground state. However they are good candidates for the (e, e'7) 
study where the L /T interferences can be isolated and the matrix elements of the 
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electric transverse multipole computed. Introducing the ratio between the electric 
transverse multipole and the Coulomb one 



^el 



(If II A || h) 



(5.4) 



and considering that £a ^ 1 5 there will be considered only the linear terms in 
the angular distribution function. The transversal-longitudinal interference term, 
which is linear in £ A is the most interesting. It can be isolated performing mea- 
surements at 7 = 0° and 7 = 180° with unpolarized electrons 



d 2 a 



dQ e dQry 



nepol 
4>-,=0° 



d 2 



a 



dVt e dVL^ 



nepol 
7 =18O°/ 



pA+->(A-2)H 
_7 

total 



^TraMottf, 



-1 



q f 2 + tan 2 °-G T Ux ((If II M x (q) \\ h)) 



where 



G 



V2 



'A + l 



TL 



7(2A - 1) V A 



5P 2 1 (cos# 7 ) -3 



' A + 2 N 
2A-3 



P^cos 7 ) 



(5.5) 



(5.6) 



Therefore, chosing convenient values which maximize (5.5), it is possible to deter- 
mine (If || M\(q) || || T x (q) || ij) for different values of g. 

Let us consider as an example the study of (e, e'j) processes when the excited 
states are vibrational collective modes. In the incompressible liquid drop model 
the density operator is given by 



Pn(t) 



3eZ 



e 



and the current density operator by 

t t \ 3eZ 1 • 
JN{r) = ^Rl^l aim 



R (l + J2 a lm Y l*m) ~ r 

Ira 



(5.7) 



V — Y lm (0A) 

K 



O(Ro-r) 



(5.8) 



For the transition to the one-surfon 2 + state, the Coulomb multipole is 

3eZ I 5 



(2+ || M 2 (q) || 0+) 



4tt v 2(B 2 c 2 y/ 2 



32(qRo) 



(5.9) 



where B 2 and C 2 are the inertia and rigid parameters. Since the liquid drop velocity 
field is postulated to be irotational, the contribution of the transverse multipole 
reads 



<2 + \\n\q) || + ) = --W^(2 + \\M 2 (q) || 



uj 2 3 



-\+\ 



(5.10) 
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with co>2 ~ 36A 1 / 2 MeV. In order to take into account the magnetization com- 
ponents of the current in eq.(5.8), a crude model for the magnetization density is 
used [[14] 

Mr) = ^Mr)L (5.11) 
and the electric transverse multipole becomes 

(2+ II *?(,) II 0+) = -ffi (l - |^) <2 + || *,(,) || 0+) (5.12) 

Introducing the deviation from the Siegert theorem (4.25) it is then possible to 
establish a connection with the interference factor (5.4) as follows [£23[] 

l2o , . (2+ II M 2 (q) || 0+) 
V 3^2 (2+ || M (g) || 0+) 

The dependence of 772 on momentum transfer q is given in Fig. 13, for 16 and 90 Zr 
in both cases : irotational flow and with nonzero magnetization ( \i = 0.5 ). Notice 
that for 16 the magnetization contributions have a sensitive effect : the deviation 
from the Siegert theorem is enhanced. For the nucleus 90 Zr, the magnetization 
contributions are less important. 

As we saw earlier the deviation from the Siegert theorem increases when we 
pass from the IF model to the RR model. Although the microscopic nature of the 
magnetization current is different from that of the RR current from geometrical 
point of view they have the same rotational structure leading to the same effect : 
enhancement of the toroidal transitions. It can be also noticed that the deviation 
from the Siegert theorem is stronger in heavy nuclei than in light nuclei. This last 
conclusion is in agreement with the statement that we have made on the electric 
transverse form factor of DTM, i.e. the toroidal effect is intensified in heavy nuclei. 



6 Conclusions 

As we mentioned in the begining the scope of this work was the study of collec- 
tive motions with a toroidal electromagnetic structure or to extract the toroidal 
contribution in nuclear transitions with a mixed rotational-vibrational spectrum. 

Considering the giant isoscalar resonance l - DTM, we calculated for a group 
of spherical nuclei ( 40 Ca, 90 Zr and 208 Pb ) the transverse form factors and we 
underlined the shift of the toroidal effect toward small momentum transfers when 
we pass to nuclei with large Z and A. Simultaneoulsy we testify an enhancement 
of the dipole toroidal response in heavy nuclei. The probability to excite the 
DTM with photon probes is small according to the calculation we have made. 
However the possibility to use electron inelastic scattering offers the promise to 
detect DTM because in such reactions it is possible to vary independently the 
momentum transfer and the excitation energy along with the scattering angle 9 
dependence of the differential cross section. Since the DTM is a transversal flow 
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associated to the Hill ring vortex the electroexcitation at 180° angles seems to be 
favourable since other modes which could occur, being especially of longitudinal 
nature, are suppressed. We have also extracted the longwavelength limit of the 
electric transverse form factor and we emphasized that it is proportional with the 
transition toroidal dipole moment. This dynamic characteristic associated to the 
DTM for small momentum transfers depends linearly on aZ, which explains the 
smallness of toroidal effects for nuclei with small numbers of protons. 

Another type of resonances that we studied are the spin dependent (spin - flip) 
modes. We emphasized their purely toroidal electromagnetic structure and their 
possible investigation with leptonic probes at 180° scattering angles and large mo- 
mentum transfer. These conclusions stems on the calculation of (e, e') differential 
cross sections of Goldhaber - Teller resonances and 1~ electric spin - flip modes 
(s-is) in 12 C, 16 0, 40 Ca and 208 Pb. 

In the second part of this work we foccused on the study of toroidal contribu- 
tions in the excitation of g.s. band of even-even nuclei from the rare-earth region 
( 152 Sm, 166 Er ). We introduced a quantity which describes the deviations from 
the Siegert theorem or in other terms the contribution of higher order terms in the 
momentum transfer of the electric transverse multipole, i.e. the toroidal multipole 
moments and their mean square radii. For the RR which is a submodel of the 
Riemann rotator the deviation from the Siegert theorem has a slope larger than 
that of the IF model. Based on this observation we showed that defining the real 
electric transverse form factor and taking the first order approximation in the mo- 
mentum transfer we can extract the transition toroidal quadrupole moment in the 
same manner in which we determine the charge or magnetic mean square radii : 
calculating the slope of the real Coulomb (magnetic) form factor. 

We computed the transition toroidal quadrupole moment and we determined 
that it depends smoothly on the vorticity This observation strengths our opin- 
ion that the toroidal moments of arbitrary multipolarity give a measure of the 
intensity of vortical electromagnetic currents in the same maner the charge multi- 
pole moments are associated to charge distribution or irotational electromagnetic 
currents. 

The electroexcitation differential cross sections of 2 + and 4 + levels from the 
g.s. band of the RR model may be approximated taking into account taking into 
account the toroidal form factors only, for a wide range of momentum transfer. This 
result emphasize the importance of toroidal multipoles relative to the Coulomb 
multipoles at 180° angles, regardless of the momentum transfer. 

We showed that in coincidence (e, e'7) reactions, by separating the longitudi- 
nal/transversal interference term we can measure the deviation from the Siegert 
theorem in the IF model and with non-zero magnetization components of the nu- 
clear current. For light nuclei the deviation from the Siegert theorem is enhanced 
by the existence of magnetic currents. 
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